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We study gravitational lensing by quadrupole potentials within the linearized gravity ap¬ 
proximation and the integration over the unperturbed photon trajectory. It is well known that 
the quadrupole potential contribution to the deviation angle is much smaller than that of the 
monopole one. We show, however, that quadrupole potentials can change the photon polar¬ 
ization vector, but there is no contribution from the monopole term to the first order. The 
effect is maximal when the axis of the quadrupole is tangential to the photon trajectory and 
it is proportional to the frequency (rate) of the quadrupole moment of the deflector. The sec¬ 
ond order correction of the monopole potential to the polarization is canceled away by the 
renormalization of the polarization vector. 


PACS: 04.25.Nx, Post-Newtonian approximation; perturbation theory; related approxima¬ 
tions 


1 Introduction and motivation 

The common belief in physics is that Einstein’s general theory of relativity inevitably implies 
the appearance of gravitational waves. The mathematical structure of these waves is defined as 
quadrupole radiative fields. 

However, one can question about two issues relevant to gravity and gravitational waves: (A) 
The existence problem: to prove the existence of gravitational waves, it is necessary to derive 
the exact wave equation from general relativity, but this is impossible because of the presence of 
at least Newtonian ’’monopole” terms that spoil the structure of the ’’wave” equation [1—4]; any 
approximate or even exact solution of Einstein field equations cannot isolate or neglect Newto¬ 
nian from quadrupole terms. (B) The locality problem: assuming the presence of gravitational 
waves infers that a decomposition of the total tensor field contains terms affecting and imply¬ 
ing the curvature of spacetime, essentially the nonlocal physical process, and at the same time 
the local radiative quadrupole field as gravitational waves; such a decomposition is difficult to 
comprehend both physically and mathematically. 
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The classical electrodynamics was the inspiration to introduce gravitational waves, but (a) 
the dipole radiative field is a solution of the exact Maxwell wave equations and (b) the electro¬ 
dynamics is completely defined only by the local gauge field. As a consequence, the predicted 
dipole radiation has been verified experimentally. 

It is possible to understand the energy loss of binary systems as the effect of quadrupole 
potentials and not as a result of radiated gravitational waves [5], 

In this paper we inspect gravitational lensing by quadrupole potentials which are parts of the 
metric tensor, such as monopole terms, affecting spacetime curvature. It should not be confused 
with the lensing of the assumed gravitational radiation by cosmological environment or some 
astrophysical sources that one can find in the literature. 

2 Lensing equations 

Gravitational lensing by quadrupoles is an old subject of investigation [6-9], but the first reliable 
calculation is due to Damour and Esposito-Farese [10]. They have shown that the incorrect 
results of Ref. [11] arise from a naive ’’plane wave” approximation. Unfortunately, the same 
unreliable approximation is used in Ref. [12]. 

Let us briefly summarize the short-wave approximation of Maxwell equations in curved 
spacetime and the transport along the rays of the wave vector (t 1 '), the scalar amplitude (a) 
and the complex polarization vector (V v ) [13]: 

t = —, £% = 0, V a = £ 0 V a . 0 =0,d+ -a£% = 0, (1) 

= diag(- 1, +1, +1, +1), U;U" = 1, i a i a = 0, 

a = v a r = 0 

We can write the transport equations more explicitly using Christoffel symbols: 

^ = \re%g af3 , (2) 

~ = ~(V + r^ A ), (3) 

a at; 2 

^ = rVpT^. ( 4 ) 

The linearized gravity approximation and the integration over the unperturbed photon trajec¬ 
tory [13] of the transport equations, lead to the following lensing equations: 


g^u = + h,j. v , b a = impact vector , = 0, 

1 f + °° 

Mp = out ) - £fj,{in) = -rt 0 / d£d^h a p(b x + ££ x ) 

^ J — OO 


( 5 ) 
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/ ,\ -j ,+oo 

Ain a = ln a ° U = — - l x / d£d\h(b M h = r] a,3 h a g , (6) 

a{in) 4 7,^ 

i r + °° 

AV a = V a (out) - V a (in) = - i> 3 V^ / dtidahfr + d 0 h ai - d^W + ^)-(7) 

" J — oo 

A direct evaluation of the lensing equations with the monopole potential straightforwardly 
gives 



2G n M 


Qcx.fi 


A In 0 = 0, = 0, 

= any spacelike vector orthogonal to 


( 8 ) 


Thus, besides the resulting nonvanishing Einstein deflection angle [13], the monopole poten¬ 
tial does not influence the transport of the scalar amplitude or of the polarization vector to the 
first order. 

For the calculus with quadrupole potentials, we adopt the formalism of Damour and Esposito- 
Farese [10] performing the integration in the Fourier space. They evaluated the lensing of the 
wave vector, which we therefore omit and concentrate on the lensing of the scalar amplitude and 
the polarization vector. 

One can read the Fourier transform of the potentials where the energy-momentum conserva¬ 
tion is also included [10]: 


h^(x x ) = 


d 4 k 


)e ' x = k v x„,T = 


hit V (k x ) = 167tGat 


(27r) 4 

% v (k x ) - \r lllv f{k x ) 


kP. 


k l k 3 


T °‘ ~ k° Tij ' T °° ~ (k 0 ) 2 ^' ^ 


k ■ k — iek° 

We can fix the constant wave (f^), impact (&") and polarization (V a ) vectors: 


r = (i,o, o,i), v* = (o, 6 ,o,o), = (o, Vi, 14,o). (io) 

There is no influence of the quadrupole as a deflector to the scalar amplitude: 

A In a(quadrupole) = 0. (11) 

Changing the integration variable k =| k | to u = \Jk 2 — to 2 [10], one obtains the expression 
for the lensing of polarization: 


k M = (u>, k sin t? cos (j>, k sin t? sin (j>, k cost?), k ■ V = k v V v , 
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j(~i r+oo r +oo p2n 

AV a =-g- dco du d<f)ue ibucos< t'(u 2 — isu ) -1 

^ J — oo J 0 </ 0 

x{fc ct Ty(-fc : 'V i a; _1 + <5j 3 Vj) - (fc • V)Ty [<5 a o(fc®/^w -2 - 5 i3 fc-'w _1 ) + + &,- 3 )] 

+^(fc • V)£ a (-k i k i w- 2 f ij +f ii )}. (12) 

Neglecting the k dependence (quadrupole approximation) of the deflector field [10] 


Tij(u,k) ~ fij(u, 0) = (13) 

Dij(u}) = J d 3 xx l x^ T°° (u, x), Dij{to) = J dte lut £>y(i), 

the rest of integrations could be performed with elementary integrals but with a careful regu¬ 
larization procedure [10]. 

The result for the lensing of polarization is of the form 


AKtnt dG N V 2 dD 12 (t) . . 4G n Vi dDi 2 (t) 

AVi(Q) = —p- 37 — \t= o, AV 2 (Q) = - ^ - 73 — |i= 0 > 


dt 


b 2 


dt 


(14) 


A V 0 (Q) = - 


2 G Nnr ,„ , n , n , *- 2 ^n(t) 

[ l 7 ! (£> 22 ( 0 ) - £>11 ( 0 ) + —- 377 - | i =0 


b 3 

14(2£i2(0) + 6 : 


2 d 2 D 12 {t) 


dt 2 


2 

U=o)], 


<9f 2 


6 2 d 2 D 22 (t) 
T dt 2 


t=o) 


AV 3 (Q) = - 


2G » [Ktft.fO) - D II( 0) - £***> lt=0 lt=0 ) 


b 3 

V 2 (2D 12 (0) - b 


2 d 2 D 12 (t) 


dt 2 


2 <9f 2 

U=o)]- 


<9f 2 


One can easily verify the following relation 


V^AV^iQ) + V^AV*(Q) = 0. 

From the exact integral for the evolution along the geodesic (C) of the polarization vector 

AV^ = [ dt(°V fj T%, (15) 

Jc 

one can easily deduce the second order corrections perturbing Christoffel symbols [2], wave and 
polarization vectors, as well as the photon trajectory with regard to the monopole potential: 


AV a ( M ) — X a + dJa + /C a + £- a + Ada, 


(16) 
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2a = \/V X d^h^{d a hf K + dphW - d K hW), 
hfy = ^(i )«7 = —rf' 1 rf‘ K h£l, Cq = unperturbed trajectory 


(17) 


Ja = 


h (2) ~ - 
'*’00 — 


\t p v« [ dt{d a h$ + dph,™ d K h&), 

^ J C n 

) h^> - n 
2 h n iO ~~ u > 


>Co 
^2 7 i ,r2 


u ( 2 ) _ G* n M\ s 

I ™ 12 > I zt 12 


(18) 


/c Q = V [ diM^[d a hf K + 

j Cq 

A4(0 = d£d a h$(b x + ^ A ), 

^ J — oo 

. . /.. „ —2M£ 2M „ 2M . 

=>■ A£ a (^) = Gtv( 0 , —-^=- 7 —, 0 , — , ), 

Ao = >/ deAvwg + a^-a^g), 

*2 Co 

a 14 = ^ de^/ig!++en, 


(19) 


( 20 ) 


AFa = Gjv(—T —(1 + + ? 2 ), ■ 7_7 7_7 

& #+F 5 


A/ i'i A/V 2 MVi 


Ma = 

Cl = 
= 

Ax a = 


[ dt(d a h% + dphW - 0 */$). 

« Ci 

perturbed trajectory , 

J dZAx^ldahW + d fi hU ~ d K h§], 

f ^{d a hf K + - 0 */$), 


(i + ^/V^T^)), 


( 21 ) 


Acc a = Gjv(M(1 + U\/b 2 + e), 


3bM 

V&re 


, 0,3M(1 + £/\A 2 + ? 2 ) - 2Mln i±xE±I!) ) 

e = small positive real cut — off. 
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The coordinate singularity in Eq.(21), regulated by the cut-off e, disappears after the inte¬ 
gration over the unperturbed photon trajectory, as one should expect for an observable. After 
performing integrations, only one term contributes to the deviation of the polarization vector: 


X Q 

Joe 

JC a 

c a 

M a 


— G 2 N (- 

= g 2 n {— 


G%( 0 , 


2 ttM 2 , 


b 2 

7 tM 2 


-Vi,0,0,- 


2ttM 2 


-Vi,0,0, 


b 2 
SttAI 2 


Vi), 


Vi), 


2 b 2 ’ ’ 8b 2 
4M 2 4M 2 2t tM 2 

Jj2 12 ^ 2 ’ ^2 


b 2 b 2 

2 3irM 2 

G At( ora 


— g 2 n {— 


2 b 2 
2 t rM 2 


26 2 


6 2 


-Vi, 0 , 0 ,— 


2irM 2 


b 2 


Vi), 


AVi(M) = AG \ 2 I2 Vu AV 2 (M) = AG i M2 V 2 . 


( 22 ) 


52 - ^ ' b 2 

However, this contribution is canceled away by the renormalization of the polarization vector: 


o /^2 M 2 

AVjx = AV^Q) + AV„(M), (V"* + + AVj,) = 1 + + 0((AH) 2 ), 

V* = (1 + + AH M ), = 1 , 

4X 12 AT 2 

= yM-A-+ AV" + OffAH) 2 ) = V" + AV^Q) + OffAH) 2 ), p = 1, 

b 2 

3 Results and discussion 

For example, take a binary system and its quadrupole moment when a surface of motion is in the 
xy plane and the z axis is defined by the direction of the incoming photon. In this case, the deflec¬ 
tion of polarization is maximal. Generally, the Euler rotation matrix R projects the quadrupole 
moment of arbitrary orientation to the frame defined by the wave vector of the deflected photon 


Dij (f) 


m 1 m 2 
m 1 + m 2 


cos 2 (f H) 2 sin( 2 f 2 f) 0 

2 sin(20f) sin 2 (Of) 0 

0 0 0 


O = 


2-7T 
~P ' 


P = orbital period , 


circular orbit : x = dcos(Of), y = dsin(Of), 

Dij(t) = R{0E,4>E)3x3£ ) ij(t). 


(23) 


The complex polarization vector and its deflections could be cast into a more conventional 
form to recognize linear, circular, or elliptical polarization states (here f3 is the phase-difference 
between two independent polarization states) [14] 
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V = Vix 0 + V 2 y 0 = V-e + + V+c. L, 

v± = ± iV 2 ), I I 2 + I 14 | 2 = 1 , e± = -L{So ± iyo), 


V+ = P+ , V-= P -e*, p = — = 


1 ~P\ 


P+ P+ 

==> 3?Vi = — =(p+ + p- cos/3), SVi = ~/=P- sin/^i 
V2 v 2 

5RV2 = —t=/ 9- sin/3, 9^2 = -4=(p- cos/3 - p+), 

V 2 v 2 

Ap+ = -i=(RAVi - 9AF 2 ), (24) 

v 2 

A/3 = -(QfAVi - 5RAV 2 + v / 2sin/3^-Ap+), (25) 

V2p-Cos/3 P- 

quadrupole contribution{Eq.(lA)) : A/3 = ^ |t=o- Ap+ = An = 0. (26) 

0“ ot 

Finally, let us present some numerical estimates of the effect assuming high-frequency (high- 
rate) deflectors: 


binary neutron star system : m = mim 2 /(mi + m 2 ), D 12 (0) = md 2 tt, 

Kepler's third law : d = fG 2 / 3 (Gjv( m i + m 2 )) 1//3 , 
to ~ M 0 , d ~ 10 _4 s, O ~ 10 3 s _1 , b ~ 10 _3 s => A/3 = 0(1CT 4 ), 

and similarly, for the accretion by supermassive black holes or coalescing black hole binaries. 
Although the effect is small, the deviation of polarization (change of the phase-difference 
between the two independent states) by the quadrupole is free of any ’’background” monopole 
contribution and it is worth measuring. 
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